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ABSTRACT 


Crack growth and closure behavior of thin and thick middle-crack tension 
specimens under constant amplitude loading were investigated using a 
three-dimensional elastic-plastic finite-element analysis of fatigue 
crack growth and closure. In the thin specimens the crack front closed 
first on the exterior (free) surface and closed last in the interior 
during the unloading portion of cyclic loading. The stabilized 
crack-opening stresses of thin (t=4.78mn) and thick (t=12.7mm) 
middle-crack tension specimens are determined at the interior and 
exterior regions under constant amplitude loading (R=0.1, Sj Te:x /a o = 0. 25, 
o i s the effective yield stress). For the thin specimen the stabilized 
stress levels were found to be 0.34 S max and 0.56 Sj—™, whereas for the 
thick specimen the corresponding stress levels were 0V26 S max and 0.46 

®max, respectively. 

A load-reduced displacement technique was used to determine crack-opening 
stresses at specified locations in the plate frcm the displacements 
calculated after the 7th cycle (using unloading and reloading portions of 
cyclic loading) . All locations were on the plate exterior surface and 
were located near the crack tip ( about 0 . 9 nm) , behind the crack tip 
(about 4.30 nm) and at the centerline of the crack. With this technique 
the opening stresses at the specified points were found to be 0.52, 0.42 
and 0.39 times the maximum applied stress. 
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NOMENCLATURE 


b ha If -width of specimen, 

c half-length of crack, 

half-length of initial crack, 
d smallest element size along crack plane, 

E modulus of elasticity, 

h half -height of specimen, 

R stress ratio (ratio of minimum to maximum applied stress) 

S applied stress, 

S Q crack-opening stress, 

Sflax maximum applied stress, 

t half-thickness of specimen, 

t^ thickness of i^ layer, 

V displacement in y-direction, 

effective yield stress, 


nm 

TTTT1 

nm 

nm 

MPa 

nm 

MPa 

MPa 

MPa 

nm 

nm 

nm 

MPa 
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1 INTRODUCTION 

The growth of a surface or through crack in a structural component is a 
three-dimensional process. Near regions vhere a crack intersects a free 
surface, the material is under plane-stress conditions. However, the 
crack front in the interior of the body nay be under plane-strain condi- 
tions. Previous experiments and analyses have shown that the crack 
closure behavior is greatly affected ty the state of stress. Closure 
under plane-stress conditions is much greater than that under 
plane-strain conditions. 

During the past two decades many researchers have used Elber's crack- 
closure concept [l-3] to explain various experimentally observed fatigue 
crack growth phenomena and also for prediction of fatigue crack growth 
under various load histories. Several researchers [4-14] have analyzed 
the behavior of crack growth and closure of specimens under various load 
histories under both plane stress and plane strain conditions using two 
dimensional elastic-plastic finite-element analysis. 

However, the crack closure behavior of finite-thickness plates under 
cyclic loading is not fully understood and two-dimens iona 1 finite-element 
techniques cannot be used for analyzing these specimens. In order to 
analyze crack growth and closure behavior of finite thickness plates, a 
three-dimensional elastic-plastic finite-element technique was developed 
[l5]. The model was used to determine crack-opening and closure stresses 
of a middle-crack tension (MT) specimen both under constant and variable 
anplitude leading [ 16 ] and also to study the influence of specimen 
thickness an the conputed crack opening stress levels [l7 ] . 

The purpose of the present study is to investigate crack growth behavior 
of thin and thick MT- specimens under constant anplitude loading 
conditions. The stabilized crack opening stresses of interior and 
exterior regions and the closure and opening profiles of the 
crack-surface plane after the 10th cycle are presented for both thin and 
thick plates. The effect of thickness on crack opening stress levels, and 
plastic zones of interior and exterior regions is demonstrated for the 
thick specimen after 10 cyclic loadings. 
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To determine the crack -opening stresses at specified points cn the thin 
plate exterior surface, a load-reduced displacement technique was used 
based cn displacements calculated after the 7th cycle (using unloading 
and reloading portions of cyclic loading). This technique is similar to 
the compliance techniques used by most experimentalists to measure crack 
opening load levels. All analyses are performed under constant amplitude- 
loading (R=0.1) assuming elastic-perfectly plastic conditions . 

2 NUMERICAL ANALYSIS 

The three-dimensional elastic-plastic finite element procedure used to 
simulate crack extension and crack closure has been described in detail 
elsewhere [l5]. Here, only a brief sunrriary is given. 

The middle-crack tension specimen (Fig.l) is modelled with linear 
eight -nod ed isoparametric solid elements (Fig.2). Due to symmetry 
conditions, only one-eighth of the specimen was analysed. 

2.1 Elastic analysis 

The small-strain theory, linear-elastic material properties and a 2x2x2 
quadrature scheme were used to calculate the element stiffness natrix 
[l8]. The equations of equilibrium were obtained frcm the principle of 
minimum total potential energy and are expressed as 

[K e ] {u} = {p} (1) 

where matrices K , U and P represent total structural stiffness, nodal 
displacements, and applied loads, respectively. Equation (1) was solved 
using Cholesky's factorization technique. The stresses and strains were 
evaluated at the eight quadrature points in the element. 

2.2 Elastic-plastic analysis 

Hie elastic-plastic analysis was based cn Zienkiewicz et al., initial 
stress method [ 19 ] . Yielding vsas evaluated at every quadrature point 
using the vcn Mises yield criterion. Stresses in excess of the yield 
stress (for an elastic perfectly plastic material) were redistributed 
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using Drucker's normality flew rule [20 ]. TOe analysis used snail-strain 
and incremental plasticity theory. The complete algorithm is given in 
ref. [l9]. The rain advantage of this method is the repeated use of the 
factorized elastic stiffness matrix. The equations of equilibrium which 
govern the elastic-plastic defomation are 

[kJ{u}={p}+{q} (2) 

Where Q is the plastic-load vector. The Q forces maintain the plastic 
strains in the elements. In the initial-stress method, the solution to an 
elastic-plastic problem is obtained by applying a series of small load 
increments to the structure until the desired applied load is reached. 

The load increment was chosen as 20% of the load required to yield the 
first element. At each load increment, eg. (2) was solved iteratively. The 
iteration process was repeated until the effective stress was within 2% 
of the yield stress in all yielded elements. Usually 5 to 15 iterations 
were required for convergence. 

2.3 Crack extension and closure analysis 

Fatigue— crack growth simulation involves opening and closing of the crack 
surfaces, Wiich leads to a variable boundary- value problem . in ref. [ 4 ], 
an algorithm is given to modify the factorized stiffness matrix to change 
the appropriate boundary conditions. The technique involves introducing a 
very stiff spring to restrain a boundary displacement car removing the 
spring to free a displacement restraint (to epen the crack). The "rigid" 
spring stiffness ( 10 ^ times the modulus of elasticity) was added to the 
diagonal coefficient in the conventional elastic stiffness matrix. To 
extend the crack, the stiffness of the springs at the crack front was set 
equal to zero in the modified stiffness iratrix K g and the forces F that 
were carried by the springs were then applied to the crack- front nodes. 

The equilibrium equations are 


[k 6 ]{u}={p}+{q}+{f}. 


(3) 
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The nodal forces F, used only during crack extension, ware arbitrarily 
chosen to be released in five equal load steps to advance the crack 
front. To insure that the stress and total strain increments in the 
adjacent elements satisfied the appropriate yield conditions and flew 
rule, the iterative analysis [ 19 ] was repeated and the forces previously 
carried by the broken springs were released and redistributed. 

During cyclic loading, the nodal displacement along the crack plane were 
monitored to determine whether the crack surfaces at any node ware to 
open (positive displacement) or close (negative displacement). If the 
crack surfaces should open, the stiffness of the boundary spring would be 
set equal to zero and the stiffness matrix would be updated. If the crack 
surfaces should close, the spring stiffness would be set to an extremely 
large value and the stiffness matrix would again be updated. No attempt 
was made in any of the analyses to include a failure criterion for crack 
extension. The crack was extended one element size (0.03nm) per cycle 
when the maximum applied stress was reached. 

3 SPECIMEN CONFIGURATION AND LOADING 

A middle-crack tension specimen (Fig. 1) was studied under constant amp- 
litude loading (R = 0.1, S /o = 0.25) in all the analyses. The 
material was assumed to be elastic-perfectly plastic. The dimensions of 
the specimen were b = 38.1 rrm, h = 2 b and t = 4.78 nm and 12.7 nm, re- 
spectively. The crack was extended by one element size (0.03 nm) at the 
maximum applied stress of each load cycle. The initial crack length c^ 
was 18.6 nm. The nodulus of elasticity E was 70 000 MPa, Poisson's ratio 
was 0.3, and the effective yield stress was 345 MPa. The finite-ele- 
ment idealization of the model is shewn in Fig. 2. Because the specimen 
configuration and loading were synmetric, only one-eicht of the plate was 
modeled. A four-layer model was used in the present analysis. The layer 
thicknesses (t^) were diosen as 0.5 t, 0.25 t, 0.15 t, and 0.1 t where 
the smallest layer was located on the exterior (free) surface of the 
specimen (z = t). The mesh had 3804 elements and 5145 nodes with 15435 
degrees-of-freedom. For purposes of discussion, the 2=0 plane is referred 
to as the "interior" plane and the Z=t plane as the "exterior" plane. 
Figure 2 also shews elements on the "interior layer" and "exterior 
layer". 
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4 RESULTS AND DISCUSSION 

First, the analysis was applied to the thinner specimen (t=4. 78rrm) . The 

stabilized opening stresses for interior and exterior regions after the 

10th cyclic loading, already reported in Ref. [l6], were found to be 0.34 

and 0.56 times the maximum applied stress (S ) , respectively, as shown 

nvnx 

in Fig. 3. Next, the same finite-element mesh was used and the thickness 
of the layers was increased. Ihe specimen thickness was 12.7 nm and the 
model was exercised under the same type of loading which was used for the 
thinner plate. The stabilized crack-opening stresses of interior and 
exterior regions after 33 cycles were found to be 0.26 and 0.46 times the 
maximum applied stress (S max ) , respectively, as shewn in Fig. 4. 

The drop in opening loads can be explained as follows. The constraint 
effect at the interior region (z=o plane), causes a variation in 
displacements through the specimen thickness which gives different 
opening loads for interior and exterior regions, respectively. 

The interior region (z=o plane) is considered to be under plane strain 
conditions, whereas the exterior region (z=t plane) is approximately 
under plane stress conditions. The crack-opening stress level of the 
interior region first reached a maximum peak (S /s =0.34) and then upon 
further crack-ex tension dropped to a ratio of 0.28. A similar trend in 
crack opening stress is reported fcy Fleck and Newman [ 14 ] under 
plane-strain conditions. The exterior region crack-opening stress level 
experiences a similar behavior and dropped to a ratio of 0.52. The 
constraint effect which gives different deformations on the interior and 
exterior regions of a three-dimensional model, aauses a variation in 
opening stresses through the specimen thickness. The crack -opening 
stresses of thin and thick specimens through specimen thickness are shown 
in Fig. 5. The thinner plate crack-opening stress level of the interior 
region is less than that of thick plate, whereas the crack-opening stress 
level of the thin plate for the exterior region is higher than that for 
the thick plate. 
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Hie closure and opening profiles of the crack-surface plane for the 10th 

load cycle of the thinner specimen, already reported in Reference [16 ], 

are shewn in Fig. 6. The interior region is fully epen at the stress 

level of 0.36 At a stress level of 0.56 the surface nodes 

behind the crack-tip node are still closed. As socn as the applied stress 

reaches beyond 0.56 S , the entire crack -plane is fully open. 

max 

The closure and opening profiles of the crack-surface plane for the 10th 

load cycle of the thicker specimen are shewn in Fig. 7. At the stress 

level of 0.34 S , the interior is fully open while the surface nodes 
max 

behind the crack-tip node are still closed. As the applied stress reaches 
beyond 0.34 S , the entire crack surface plane is fully epen. Clearly, 
the thicker specimen has a smaller surface layer influenced by the state 
of "plane stress" at the free boundary, due to more plastic constraint, 
and consequently the overall crack closure behaviour is more dominated ty 
the interior plane strain region. In this analysis, the thick specimen 
had the same thickness ratio (thickness layer over plate thickness) as 
the thin plate. To obtain a better result, the thick plate with several 
layers should be used in the analysis. 

The effect of constraint can also be visualized ty comparing the size of 
the plastic zones resulting from different deformations through the 
specimen thickness after the 10th cyclic crack extension for the interior 
and exterior regions of the thick specimen as shewn in figs. 8 and 9, 
respectively. 

4.1 Load-reduced -displacement compliance technique 

The load-reduced-displacement technique is frequently used ty experi- 
mentalists to obtain the crack-opening stress from displacements measured 
at various locations along the crack plane. Here, the load— reduced- 
displacement technique was used to derive crack— opening stress levels 
from displacements at specified locations A^, A 2 and A^ along the crack 
plane on the outer surface far the thin specimen ( t=4 . 78 rrm) , (see 
Fig. 10). The Aj^ location was first chosen ty Elber [ 3 ] to make his 
experimental measurements. Location A 2 is 4.30 nm behind the crack front. 
The A 3 location is at the center of the specimen which is conntxily chosen 
for experimental measurement of crack -opening stress levels. 
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The applied stress is plotted against the crack-opening displacement for 
the specified points A^, and A^ as shewn in Fig. 10. These curves are 
constructed using the displacements calculated for 65 load increments. 
Points A 2 and A^, Which lie on the crack surface, give larger 
crack-opening displacements than point Al. However, the displacements at 
point should be more sensitive to the closure behavior on the outer 
region of the model. Figs. 11 to 13 shew the load-reduced displacement 
plots for points A^, A 2 and A^, respectively. These load-reduced 
displacement plots were obtained using the following steps. 

A straight line is fitted through unloading points of the 7th cyclic 
loading. Knowing the equation of this line, an expression for the 
displacement quantities is obtained. For each load increment, the 
corresponding displacement quantity is then found from the equation of 
the straight line. The difference between the finite element calculation 
and the straight-line displacement is called the reduced displacement 
which is due to the effects of plasticity. 

The applied stress and reduced displacement quantities were found using 
the unloading and reloading portions of the cyclic load history. For 
unloading and the range of applied stress between 86.2 MPa and 68.95 MPa, 
the reduced displacement values are zero. Between 68.95 MPa and 41.37 MPa 
stress range, the reversed plasticity at the crack front causes the re- 
duced displacement values to decrease. Below 41.37 MPa, reversed 
plasticity and crack closure causes erratic behaviour of the reduced 
displacement plot. Crack closure tends to increase the reduced 
displacement value (opposite to that caused by reverse plasticity). At 
the minimum applied stress, S . =8.62 MPa, the naximum reduced 

displacement value was obtained. For the reloading portion of the load 
history, the curve goes to the left because the crack-surfaces are 
opening behind the crack front. For the point at which ®S/ 3AV = the 
crack-opening stress was obtained. The crack-opening stress for point A^ 
was found to be 45 MPa. Upon further reloading, the curve moved to the 
right and came back to its starting point. 
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A similar trend was found for points and A^ . The crack -opening 
stresses for points A ^ and A^, using the load-reduced displacement plot, 
were found to be 36 and 34 MPa in Figs 12 and 13, respectively. 

The normalized crack-opening stress, as a function of normalized distance 
from the center of the crack for points A^ , and A^ is shown in Fig. 

14. Based cn these results, it can be seen that point A^ has a reduced 
displacement crack-opening stress close to the outer region reduced 
displacement crack -opening stress. Cn the other hand, point A^ has a 
reduced displacement crack-opening stress level closer to the inner 
region value . The reduced displacement crack -opening stress for point A^ 
is located between the crack-opening stresses of points A^ and A^. 

5 CONCLUSIONS 

The three-dimensional model was applied to middle crack tension specimens 
with two different thicknesses to investigate crack growth and closure 
behavior under constant amplitude loading. A load-reduced displacement 
technique was used to calculate the crack opening stresses at specified 
locations on the crack surface plane. Results frcm the analysis support 
the following conclusions : 

i) The crack opening stress level in the interior region is lower 

than that for the exterior region of the three-dimensional body 

which is due to the constraint effect. The effect of constraint on 

thickness gives a lower stabilized crack -open ing stress level of 

0.26 S and 0.46 S for the interior and exterior regions of 
max max 

the thick specimen than for the thin plate, where the same values 

became 0.28 S and 0.52 S , respectively, 
max max 

ii) The crack-opening stresses at three specified locations on the 
crack-surface plane of the thin specimen (t = 4.78 nm) were 
calculated using a load-reduced-displacement technique. The 
normalized crack -opening stress, as a function of normalized 
distance fran the center of the crack for points A^, A^ and A^ is 
shown in Fig. 14. 
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iii) Closure values determined by the load-reduced-displacement 

technique (similar to experimental studies) approach those derived 
from the point of contact at the first crack tip node in the 
finite element solution. However, due to the smaller sensitivity 
of the load-reduced- displacement technique those values are 
dependent cn location on the specimen surface. Close to the crack 
tip the value approaches the "plane stress" value, while far away 
from the crack tip a value close to the inner region "plane 
strain" result is obtained. 
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Fig. 2 Finite-element idealization of middle-crack tension specimen 
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Fig. 3 Normalized crack-opening stresses on interior and exterior planes 
under constant-amplitude loading (t=4.78 mm) 
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Fig. 4 Normalized crack-opening and closure stresses on interior and 
exterior planes under constant amplitude loading (t = 12.7 mm) 



Normalized crack-opening stress variation 
through the specimen thickness 
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Fig. 10 Applied stress against 
specified points along 
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11 Determination of crack-opening stress using a load-reduced 
displacement technique for point A ( t =4 • 78 mm) 
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Fig. 12 Determination of crack-opening stress using a load-reduced 
displacement technique for point A (t=4.78 mm) 
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Fig. 13 Determination of crack-opening stress using a load-reduced- 
displacement technique for point (t=4.78 mm) 
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